Entanglement frustration in multimode Gaussian states 
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Bipartite entanglement between two parties of a composite quantum system can be quantified in terms of tlie 
purity of one party and tliere always exists a pure state of the total system that maximizes it (and minimizes 
purity). When many different bipartitions are considered, the requirement that purity be minimal for all biparti- 
tions gives rise to the phenomenon of entanglement frustration. This feature, observed in quantum systems with 
both discrete and continuous variables, can be studied by means of a suitable cost function whose minimizers 
are the maximally multipartite-entangled states (MMES). In this paper we extend the analysis of multipartite 
entanglement frustration of Gaussian states in multimode bosonic systems. We derive bounds on the frustration, 
under the constraint of finite mean energy, in the low and high energy limit. 



I. INTRODUCTION 

The development of quantum information theory and tech- 
nologies has stimulated and motivated the scientific efforts to- 
ward the full characterization of the geometry of the set of 
states of quantum systems For both fundamental and 

technological reasons, the interest has been focused on the 
characterization of entangled states, that is, those states of a 
composite system exhibiting nonclassical correlations among 
their parties. While the case of bipartite systems has been ex- 
tensively studied and — at least for pure states — is nowadays 
well understood jst), the case of genuinely multipartite entan- 
glement is still not fully mastered. 

The present contribution deals with the problem of 
multipartite-entanglement characterization by focusing on the 
phenomenonof entanglement frustration @]. In particular, we 
will explore Gaussian states of continuous variable (CV) sys- 
tems, that is, systems of (quasi-free) quantum harmonic oscil- 
lators 101. Specifically, we consider a suitable cost function, 
the potential of multipartite entanglement, introduced in 
and extended to the Gaussian framework in i^, as a quantifier 
of frustration. The fact that this cost function cannot saturate 
its minimum value is a symptom of a sort of frustration of en- 
tanglement, induced by the geometry of the quantum phase 
space, which prevents the states to be maximally bipartite- 
entangled among all the possible system bipartitions. Here 
we consider a family of entanglement cost-functions, general- 
izing the one introduced in 121 . We hence derive new results 
on their minima in the low and high energy limits. 



II. GAUSSIAN ENTANGLEMENT 

Our analysis focuses on a system of n quantum har- 
monic oscillators, namely a set of n bosonic modes, de- 
scribed by the canonical variables X = {Xi , X2, ■ ■ ■ X2n) '■= 
{qi,q2, ■ ■ ■qn,Pi,P2, ■ ■ -Pn)- For the sake of simplicity, we 
assume all oscillators to be identical, although distinguish- 
able, with unit frequency, and set /i = 1. We follow @1 ™d 
consider the manifold of Gaussian states of the n-mode sys- 



tem. Let us recall that a Gaussian state p is characterized by 
the first and second moments of the canonical variables, that 
is, the mean (X) := tr(pX), and the covariance matrix (CM) 
V, with elements Ya.b = ^{XaXb + XtXa). We assume, 
without loss of generality, (X) = 0, and restrict our attention 
to pure states. The CM of a pure state can be written in the 
form 
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with 



T = 



X Y 
-Y X 



(1) 



(2) 



where IK is diagonal and nonsingular, O denotes the null ma- 
trix, and R is a symplectic orthogonal matrix, characterized by 
the property that U = X + «Y is unitary. Finally, we impose a 
bound on the mean energy per mode, that is, 

Vfc,fc+V„+fc,„+fc 1 \fk = l,...n (3) 
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with N the number of mean excitations per mode. This is one 
of a number of physical constraints that must be imposed on 
the system in order to make the problem mathematically (and 
physically) well posed. A different approach to entanglement 
frustration, that makes no use of energy constraints, has been 
proposed in 1.1011 . 

We consider purity as an estimator of bipartite entangle- 
ment mistl. Given a bipartition into two subsystems {A, A}, 
associated with two subsets of respectively \A\ = ua and 



1^1 



bosonic modes (ua 



n), the purity of the 



subsystems reads 



TTA = T^A 



(1/2)"^ 



(4) 



where is the sub-matrix of the CM identified by the indices 
belonging to subsystem A. Here we have assumed without 
loss of generality ^ n^. It can be easily proven that under 
the constraint ^ the minimum value of the purity is 



(1/2)'^ 



(A^+ 1/2)'^ 



(5) 
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The range of the purity is tta € [^n^^i-^)' 1]' where the value 
1 characterizes/acfo nzec/ states, and the value 7r"""(A^) char- 
acterizes those states which are maximally bipartite-entangled 
across the bipartition considered. 



A. Entanglement frustration in multimode systems 

In order to study multipartite entanglement we introduce 
the normalized potential of multipartite entanglement, a cost 
function defined for any pair (n, ua) by the (normalized) ex- 
pectation value E of the purity over all possible bipartition of 
given size: 



X{r. 



n 



E 

|A|=ri, 



(iV + 1/2)"'*E [dct(VA)"^/^ 



(6) 



The range of the cost function is contained in the inter- 
val [1, l/7r"""(^)], where the lower bound characterizes the 
so-called perfect MMES (Maximally Multipartite-Entangled 
States) I9I1, which are maximally bipartite-entangled across all 
bipartitions of size ua- However, as shown in the geom- 
etry of the manifold of CV Gaussian states prevents the min- 
imum of the cost function to saturate the lower bound, that 
is. 



mm 



[N) 



{N) > 1, 



(7) 



a feature that is interpreted as frustration of entanglement. 

Numerical results for the case of balanced bipartitions 
{ua = were presented in j^. Figures [T] |2] shows the 

behavior of X™'"^) a function of N for 7ia = [n/2], and 
for Ua — 2 < [ri/2]. appears to be a mono- 

tonic function of N. Two limiting regimes are identified, cor- 
responding to iV ^ 1 and N ^ 1: in the former case we 
notice a linear regime for increasing values of N; in the latter, 
^(n tia) saturates to a constant that depends on the value of n. 
Moreover, xfn'nA) ^PP^^rs to be a decreasing function of 
n for a given ua < [n/2], while for ua = [n/2] it increases 
with n, although oscillating between even and odd n. 

More generally, frustration — which is naturally quantified 
by X("i"i^) (^) — decreases with n at fixed ua, and increases 
with n if UA scales linearly with n (e.g. for balanced bipar- 
titions). This behavior of frustration in CV Gaussian states is 
analogous to that observed in discrete-variable quantum sys- 
tems Wl. 




FIG. 1: Minimum of the normalized potential of multipartite en- 
tanglement X(^"|n/2]) number of excitations per mode A^, for 
n = 4,5, 6, 7, 8, 9. In the region A*' ^ 1 the minimum value is 
linear in A'^, while it reaches a plateau for N ^ 1. See Sec.lIIIIfor 
details. 




FIG. 2: Minimum of the normalized potential of multipartite entan- 
glement in the case of unbalanced bipartitions xln"^!) versus A'^, for 
71 — 4, 5, 6, 7, 8, 9. In the region N <^ 1 the minimum value is 
linear in A^, while it reaches a plateau for AT S> 1. See Sec.lIIIIfor 
details. 



all modes), whose CM reads 



-YX ' YY ' 



-2r 



YYT YX^ 
XYT XX'f 



(8) 



or, equivalently, 

cosh(2r)^ 
V = 

sinh (2r) 



XXT - YYT -XYT - YXT 
-YXT - XYT YYT - XXT 



(9) 



The energy constraints read 



in. BOUNDS ON ENTANGLEMENT FRUSTRATION IN 
GAUSSIAN STATES 



In order to estimate the bounds on X^'Ji^) (-^)' we restrict 
our attention to states such that K = e''I„ (same squeezing for 



cosh (2r) < 2iV + 1 . (10) 

In the following we denote respectively X{n,nA) '^^^ X™" ^) 
the normalized potential of multipartite entanglement and its 
minimum evaluated for this family of states. 
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TABLE I: Numerical estimates of the slopes Q(„.„4) in Eq. i\6l . 
<3(n,n4) in Eq. l ll5t . of the minimum of the normalized potential 
of multipartite entanglement for N = 10, X(n','[,i/2]) = 10) 
[which approximates limjv_>oo X(n"[n/2]) (-^)]' ™d of the upper 
bound /3(,i,n_4) in Eq. il9l . The results are obtained for the case 
of balanced bipartitions, ua = [n/2]. 



n 


Q^(n,[n./2]) 


3(n,[n/2]) XT,l 


(./21,(^ = 


10) /3(n,[n/2]) 


4 


1.33 


1.333333 


1.663650 


1.666667 


5 


1.00 


1.000000 


1.332326 


1.333333 


6 


2.40 


2.400000 


2.780918 


2.795085 


7 


2.00 


2.000000 


2.203228 


2.213586 


8 


3.43 


3.428571 


5.980689 


6.074700 


9 


3.00 


3.000000 


3.470522 


3.491497 



A. The linear regime, -C 1 
From Eq. (|9]l we get 

X{n,nA) ~ 

{2N + 1)"^E [det (cosh (2r)I„^ + siiih (2r)Z^)"^/^ 



where we defined 

Z ;= 



-YXT - XY^ YYT - XX^ 



(11) 



and denotes the sub-matrix corresponding to the subset of 
modes A. Then, 



X{n,nA) 



f 2N + 1 



> E 



\cosh (2r) 

dot (I„^ + tanh (2r)ZA)"^/^ 



det (I„^ + tanh (2r)ZA) 

(12) 



where the inequality follows from the energy constraint, and 
it is saturated when cosli2r = 2N + 1. 

In the limit <C 1, we use the second-order expansion of 
the determinant, 

, N tr(M)2 - tr(M2) , „ 
det(I + eM) = l + etr(M) + — ^ ' ^ — ^ l+0{e^). 

(13) 

By setting cosh2r — 2N + 1, and noticing that tr(Z^) — 0, 
we get 

E det (I„^ + tanh (2r)ZA)"^^^ = 

1 + NE [tr(Z^)] + 0{N^^^) . 

Finally we obtain the upper bound: 

xTZa)W ^ xTCnA)(.N) 1 + TVmin {E [tr(Zi)] } . 

(14) 

It is worth noticing that the evaluation of the upper bound still 
requires a constrained minimization, which is now indepen- 
dent of N. In conclusion, in the region iV <C 1, the minimum 



X(n n 4) (-^) bounded from above by a linear function of TV. 
The value of the slope 



'Hn.riA) ■= min {E [tr(Z^)] } 



(15) 



has been evaluated numerically, for several values of n, and is 
presented in Table |T] A comparison with the numerical esti- 
mation of 



(16) 



also reported in Table U leads to conclude that the bound is 
tight. 



B. Saturation, iV > 1 

Let us rewrite Eq. (O as 



=2r 



(W + e-'^'^W') , 



where 



/ XXT -XYT 
\^ -YXT YYT 

We thus obtain 

(iV+ 1/2)1-41 



W' := 



X{n,nA) 



(e2'72)l^l 
(Af+ 1/2)1^1 



E 



YY^ YX^ 
XYT XXT 



dot {Wa + c-^^'W'a) 



(2iV + 1 - ie-2'-) 



det 



(17) 



where Wa, are respectively sub-matrices of W, W. No- 
tice that inequality (Tl\ is saturated if cosh 2r = 2N + 1. 
In the limit iV > 1 (i.e., e""" < 1), we get 



X{n,nA) 



> 2""'*E 



dct(WA) 



-1/2 



+ 0(e-2'^). (18) 



Then, by setting cosh 2r = 2iV+l, we obtain the upper bound 

min / ^ -mill 
X(ri,nA)'y^^ ) ^ X(n,nA) 

^2~"-''min|E dct(WA)"^/^ } 

=: l3(n,nA) ■ (19) 

Notice that, also in this case, the evaluation of the right-hand 
side of this inequality requires a constrained minimization, 
now independent of the energy parameter A^. These inequal- 
ities imply that the minimum X(^"Ji^)(TV) is bounded from 
above in the iV ^ 00 limit. The upper bound in Eq. (T% 
can be evaluated numerically. Table |T] shows a comparison 
between the values of X™'n4)(^)' calculated for = 10 
(where the saturation regime has been reached of all values of 
n considered), and the numerically estimates of /?(„,n^). Nu- 
merical evidence suggests that the upper bound is approached 
in the N ^ 00 limit. Comparison with Eq. (Q yields then a 
concrete estimate for the amount of frustration in the system. 



4 



IV. CONCLUSIONS 

We have presented some results on entanglement frustra- 
tion in multimode Gaussian states, quantified by the minimum 
value of the normalized potential of multipartite entangle- 
ment. Entanglement frustration arises from the impossibility 
for a multimode Gaussian state of being maximally bipartite- 
entangled across all possible bipartitions of the system. It has 
been proven in that entanglement frustration appears in 
multimode Gaussian states if n 4 ^ 2, while for qubits — for 
the case of balanced bipartion [ua = ["-/S]) — it appears for 
n = 4, 71 ^ 8 (the case n = 7 is still under debate IslfTTIl'). 

The results obtained in this note extend the numerical anal- 
ysis presented in j9^, and put it on a more solid basis, due to 
the semi-analytical calculation of the upper bounds in the low 
and high energy regimes. Our numerical analysis suggests 
that these bounds are tight. In particular, the calculation of 
a finite upper bound in the high energy regime demonstrates 
that entanglement frustration remains finite in the N 00 
limit. 

The results of the numerical analysis show a certain regular- 
ity in the estimates of the parameter a(„^„4). The numerical 
estimates reported in TableUsuggest to conjecture the follow- 
ing relations: 



a{2nA + l,nA) =nA- I 

2nA 



2nA - 1 



a(2nA + l,nA)- 



(20) 
(21) 



We remark that the upper bounds are obtained by an optimiza- 
tion over the matrices M in Eq. which have the property 
of being both symplectic and orthogonal, and hence define 
a representation of the unitary group U{n) (see, e.g., |fl2|]). 
These observations suggest that the employment of group- 
theoretical methods could lead to deeper insight into the phe- 
nomenon of entanglement frustration in multimode Gaussian 
states. 
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